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Abstract 

We argue that a strongly first order electroweak phase transition is natural in the presence 
of strong symmetry-breaking interactions, such as technicolor. We demonstrate this using an 
effective linear scalar theory of the symmetry-breaking sector. 
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1 Introduction 

The nature of the finite temperature electroweak phase transition remains an important and elusive 
problem in particle physics. The possibility of electroweak scale baryogenesis, which has helped 
generate new interest in this subject Q, requires the transition to be strongly first order; only then 
does the universe depart from thermal equilibrium, satisfying the third of Sakharov’s requirements 
for baryogenesis and avoiding the washing out of any generated baryon asymmetry. This paper 
investigates the electroweak phase transition when the symmetry-breaking sector involves strong 
interactions, such as technicolor. We use an effective scalar theory to describe these interactions, 
and find that the transition is first order if this sector exhibits a U{N)l x U{N)^ global symmetry. 
When the interactions are relatively strong, the transition is strongly first order: the discontinuity in 
the order parameter at the critical temperature is of the order of the zero temperature expectation 
value V. 
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To set the stage, we briefly review what is known about the electroweak phase transition in 
the minimal standard model. A variety of studies of the finite temperature effective potential 
have addressed this problem, using methods including the e expansion, Wilson-effective-action and 
“average action” renormalization group techniques, or graphical summation schemes with names 
such as daisy and superdaisy |^, |3| . They allow improved treatments of the infrared problems 0 
that can arise near phase transitions. The transition appears to be first order when the Higgs boson 
is lighter than the W and Z, but at most rather weakly so for an experimentally acceptable Higgs 
boson mass. As that mass increases, corresponding to stronger scalar coupling, the strength of the 
transition decreases ||^ . When the Higgs boson is much heavier than the W and Z, it is reasonable 
to neglect the gauge and Yukawa interactions when studying the transition. In that approximation 
the minimal standard model is the 0(4) linear sigma model, and the chiral transition in 0{N) 
models seems to be second order or at most weakly first order |^, |^ § . 

Beyond the minimal standard model, less is known about the nature of the electroweak phase 
transition, especially if the symmetry-breaking sector is strongly interacting, as in a technicolor 
theory. Lattice studies of QCD-like gauge theories indicate that with two light fermions, corre¬ 
sponding to an SU{2)l x SU{2)ji ~ 0(4) global symmetry, the transition is second order, whereas 
with more light fermions and a larger global symmetry, the transition is first order. This suggests 
that technicolor theories, which often do have more than two technifermions, might lead to a first 
order transition. To complement the lattice results, we analyze the strongly coupled symmetry¬ 
breaking problem using an effective scalar theory, which is partly amenable to analytic treatment 
and avoids some of the uncertainties associated with lattice techniques. 

2 The Linear Sigma Model 

The properties of the phase transition depend on the dynamics of the electroweak order parameter. 
We assign it to transform under a global U {N) lxU chiral symmetry, motivated by technicolor 
theories, in which the order parameter is a condensate of fermion bilinears. A one family technicolor 
model for example corresponds to = 8 flavors, large enough to suggest utility of the large N 
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expansion in the analysis. The most economical effective scalar theory of electroweak interactions 
is of course the nonlinear chiral Lagrangian. To allow for the existence of a phase transition 
above which the order parameter vanishes, however, we must go beyond the nonlinear theory |T^]. 
Accordingly, we will use an effective linear sigma model to describe the symmetry-breaking sector 
00- Neglecting standard model gauge and fermion interactions in comparison to this sector’s 
strong interactions, we take the Lagrangian to be 

£ = TV (Ets) - ^ (Tr Stj])' - ^TrOj . (1) 

Here we have defined the traceless operator O 2 = — {X/N) Tr (S^S), where 

£ = {(Ta + iT^a) T°' (a = 0,1,, N‘^ — 1) , (2) 

with generator matrix normalization 2 Tr = 5“^ and = X/^/2N. Higher-dimensional 

(nonrenormalizable) interactions could also be included in the effective theory, but we argue in 
section Q that their omission does not qualitatively affect our conclusions. 

In order that the tree level potential be bounded below we must have Ai > 0, and either A 2 > 0, 
or A 2 < 0 with Ai > [N — 1)|A2|. Spontaneous symmetry-breaking then occurs for negative /r^. We 
do not consider the A 2 < 0 case, in which U{N)l x U{N)r breaks to [U{N — 1) x U{1)]l x [U{N — 
1) X [/(l)]/j. For A 2 > 0 the breaking pattern is U{N)l x U{N){i — > U{N)v, and the tree level 

potential at zero temperature is minimized by a vacuum expectation value (n|E|II) that can be 

taken, after a U{N)l x U{N)ji transformation, to be real and proportional to the identity matrix; 

M 

{n\aa\n)= 0 {a = l,...,N^ -1) , 

{n\TTa\n)= 0 {a = 0,...,N^ -1) . (3) 

The A 2 independence of (n|(To|n) follows from the definition of O 2 . The W mass is given by 
so that (n I (To I If) = V = 250 GeV. The 2iV^ degrees of freedom described by the S fields correspond 
to N‘^ Nambu-Goldstone scalars vTa, and N'^ scalars (cra,cro) that describe the massive fluctuations 


3 




of the order parameter. Their masses 


2 n 2 ^2 2 


Ai + A 


2iV2 


2 

^<^0 - 2 iV 2 


2 2 
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(4) 


scale as 1/A^ relative to = 250 GeV. 

Among the Nambu-Goldstone scalars, three become the longitudinal and Z, the tto 
gains mass via anomalous breaking of the axial U{1), and the remaining — 4) are pseudo- 
Nambu-Goldstone bosons, all of which gain small masses from neglected standard model and other 
(Extended Technicolor) interactions. We neglect these pseudo-Nambu-Goldstone masses. The 
anomaly-generated ttq mass could also be implemented in our effective Lagrangian, by adding a 
nonrenormalizable (for W > 4) determinantal interaction [pi] , |^, 131, but omission of this term 
affects the masslessness of only one mode and will not change our conclusions. In the limit A 2 = 0, 
the symmetry of the theory increases from U{N) x U{N) to 0{2N‘^) and the spontaneous breaking 
produces (2Ai^ — 1) Nambu-Goldstone bosons, with only the do massive. 


The factors of N in Eq. (||) allow a nontrivial large N expansion holding Ai and A 2 fixed. In 
the A 2 = 0 limit, the leading 1/N approximation is tractable, corresponding to the familiar linear 
bubble sum. With A 2 nonzero, on the other hand, all planar diagrams involving the A 2 interaction 
contribute at leading order. Eor both interactions, the strong coupling limit sets in at Aj/IOtt^ ~ 1, 
when higher loops are as large as lower order contributions. In this limit the scalar masses of 
Eq. (^) saturate the bound 47ru/A^ Q. 


Symmetry-breaking aspects of the U{N) x U{N) linear sigma model at zero temperature in 
various spatial dimensions have previously been examined using the e expansion jllj combined with 


the effective potential |15|, and with lattice methods |16|. These investigations found Goleman- 
Weinberg [17] behavior: broken symmetry when is tuned to zero. The e expansion analysis, using 
renormalization group flow in d = 4 — e spatial dimensions, shows that a second order transition 
is not self-consistent for N > 2 [llj. However, studying the transition as a function of Lagrangian 
parameters at T = 0 is not the same as studying it as a function of temperature. Since effective 
three-dimensional theories apply to finite temperature only by assuming decoupling of nonstatic 
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modes (discussed below), they are limited to the second order or very weakly first order cases. The 


same is true of analyses of the transition strength using renormalization group flow “time” |]^]. 
Furthermore, the e expansion may break down before e = 1 or miss IR fixed points. Nevertheless, 
the above work does hint that a full finite temperature calculation might find a first order phase 
transition. 


3 Computing the Effective Potential 


We analyze the phase transition with the finite temperature effective potential defined in Euclidean 


periodic time, which describes the system in thermal equilibrium [19|. Matsubara frequencies 

= (27rnT), where n = 0,±1,±2,..., then appear in propagators. If the relevant momentum 
scales were much smaller than 2tiT (e.g. the high temperature limit), all but the n = 0 modes 
would decouple, leaving an effective three-dimensional theory. This is not true at a strongly first 
order phase transition, where we shall see that relevant scales are of order 2ttT. (The transition 
temperature will be of order v/N, while in the strongly coupled case the masses in Eq. (|^ approach 
4:Ttv/N.) We will therefore retain the full four-dimensional dynamics; we will also avoid the relative 
simplicity of the high temperature expansion, relying below on numerical integration for explicit 
results. 

We can obtain the finite temperature effective potential as the sum of one-particle-irreducible 
vacuum graphs in a background field [p^. For the linear sigma model of interest here, the calculation 


is conveniently implemented using the auxiliary field method \2A, to eliminate the Ai interaction 
in favor of a nonpropagating dimension-two field y. The auxiliary field facilitates resummation to 
all orders in Ai, at each order in 1/iV. It consists of adding to the Lagrangian a perfect square, 
which yields an irrelevant constant factor upon path integration over the auxiliary field: 

4iV2 / 


C 




(5) 


Before presenting results for the full theory, it is useful to review the nature of the phase tran¬ 
sition in the A 2 = 0 theory {0{2N'^) symmetry). Analyses using the 1/N expansion show that the 
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high temperature phase transition here is second order |^, ^]. The leading order computation 

of the effective potential corresponds to the “superdaisy” approximation |l^, | 2 ^, which can 

therefore be justified in the context of the 1/N expansion. Without superdaisy resummation, on 
the other hand, a one-loop computation using the Lagrangian of Eq. with A 2 = 0 gives a first 
order result 1^, 261. This is unreliable, since the theory’s infrared behavior is controlled by the 
effective loop expansion parameter, of order (Ai/lOvr^) T/m{a). Here rin?{a) = + (Ai/4)cr^ is 

the effective mass of the lightest excitations, in the quantum state jT) with amplitude peaked at 
the classical background field (T|(To|'I') = Na. This mass vanishes for Nambu-Goldstone bosons 
at a symmetry-breaking minimum of the effective potential, so the perturbative expansion breaks 
down due to infrared divergences, in some a neighborhood of the minimum. The 1/N expansion 
avoids this perturbative problem by a summation of graphs at each order in 1/N. 


To leading order in 1/iV, the 0(2iV^) computation involves only a single loop of the 2iV^ — 1 
Nambu-Goldstone bosons in the presence of background y and a fields, using the Lagrangian of 
Eq. dD- Because the y field is nondynamical, it can be eliminated in favor of a using the saddle- 
point equation of motion. This procedure reveals the second order behavior. It is equivalent to 
using the Lagrangian of Eq. and solving a momentum-independent Schwinger-Dyson equation 
for the Nambu-Goldstone boson effective mass, whose solution is inserted into the single loop of 
Nambu-Goldstone bosons. 


The effective potential computed in this way becomes complex for a between the symmetry¬ 
breaking minimum (where the Nambu-Goldstone boson mass vanishes) and the origin. This occurs 
for A 2 7 ^ 0 as well; in either case, the imaginary part may be interpreted |27] as the transition 
amplitude for leaving an unstable configuration, into which a classical external source cannot force 
the system. One may define a real quantity by a Maxwell construction: the “convex effective 
potential” (the convex envelope of the real part of the above effective potential) [^]. Although 
the convex effective potential approximately describes the lowest energy quantum state for a given 
(T, if a local minimum remains at the origin then the system is likely to supercool, remaining in 
the quantum state peaked at a = 0 even below the temperature Tc where the symmetry-breaking 
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minimum falls to a lower energy than the symmetric one [^|. The phase transition is thus controlled 


by the appearance and location of minima in the nonconvex potential as computed here. 


Turning now to the general case with A 2 7 ^ 0, the computation is complicated by the fact that 
planar graphs to all orders in A 2 contribute at leading order in the 1/iV expansion. Auxiliary fields 
could also be introduced to replace the A 2 interaction, but not very usefully since an infinite class 
of graphs would still contribute at each order in \/N. We instead use the Lagrangian with only 
the auxiliary y field, Eq. (I), and compute the effective potential as we did in the A 2 = 0 case: one 
loop of {2N'^ — 1) quanta in background x a fields (Fig. |T]), after which the auxiliary y field is 
eliminated using its equation of motion. This approximation sums all orders in Ai to leading order 
in 1/N. It also includes the (tree-level) contribution of the A 2 interaction to the effective masses 
of the ( 2 iV^ — 1 ) quanta in the loop, lifting the degeneracy between the iV^ Nambu-Goldstone 
modes and the (iV^ — 1 ) aa modes which are Nambu-Goldstone bosons only in the A 2 —> 0 limit. 
This approximation to the U{N) x U{N) theory should provide a reliable qualitative guide to the 
nature of the phase transition, even without contributions like those of Fig. 



Figure 1: Feynman graphs included in our computation of I4ff. The vacuum loops of {N"^) flavors 
of TTa and {N'^ — 1 ) of aa propagate in background do and x fields. 



Figure 2 : Some neglected graphs. One-loop graphs subleading in 1/A^; and higher order graphs 
suppressed by powers of 1/A^ or A 2 or both. Dashed lines are diagonalized x held propagators, solid 
lines are diagonalized a and vr fields. 
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The computation requires renormalization counterterms which we define conventionally at zero 
temperature. Diagonalizing the propagators then gives the one-loop result for the finite temperature 
effective potential as a function of a = (T|(To/A^|'I') and y = (TlxlT): 


N-^Via,x) = 




Ai 


+ \ +X)+ + 

— (T = 0 counterterms) 


( 6 ) 


where = (27rnT) and Jj, = T /(27r)“^(i^/c X)n='^oo- The first integral corresponds to the loop of 
TTa (Nambu-Goldstone) scalars whose effective mass is y; the second integral corresponds to the 
loop of aa scalars whose effective mass is x -|- (A2/2 )cj^. The subtracted counterterms are 


In -I- 


/•X+^Tr(StS) 

1 fX^ + 

X + ^Tr (StS)' 

0 

4)0 

(P + M2)2 





f X+ 

Jo 


1 fx^ + ix + ^(T^? 

4 Jo (A:2 + M2)2 ’ 


(7) 


where the subscript 0 on the integrals refers to zero temperature. The first term is an unobservable 
vacuum energy and the second corresponds to renormalization of In the final term, which 
corresponds to and Ai renormalization, we have inserted a renormalization scale in the form of 
an infrared cutoff on the integral. Physical results are not sensitive to details of the renormalization 
scheme. 

The renormalization of Ai to leading order in 1/N includes all Ai-dependent corrections, and a 
subset of the A 2 -dependent corrections. There is no A 2 renormalization in our approximation (with 
a purely no background). The corrections to Ai lead to a Landau pole at a sufficiently high scale, 
reflected in a corresponding pathology in the effective potential at a sufficiently high value of a. 
We avoid this problem by taking Ai to be somewhat less than the strong coupling limit at the scale 
Mr, which in turn we take to lie just above the momentum scales relevant in the effective potential. 
Thus the Landau pole is safely out of range at the scales of interest. 



Carrying out the summation in Eq. (^) gives 


N-^V{a,x) = 
+ 

+ 

+ 


1 2 , 

2^” ^ — 

2^4 

27r^ 

{x + ^? L(x + ^) 


X A X 3 
647r2 V 2 


J x^dxln^l —exp — + y/T^ ^ 


647r^ 
i 

27r2 JO 


M 2 


3 


2^4 /-oo / 

J X dx In I 1 — exp 


-VX 2 + 


(x + ^ 


2^2 


The auxiliary field is eliminated by solving its equation of motion 


/ dV{a, x) 

V dx 


= 0 . 


( 8 ) 


(9) 


This transcendental equation, which we solve numerically, is equivalent to a single Schwinger- 
Dyson equation for the scalar masses. Had we introduced another auxiliary field to resum the A 2 
interaction, we would have coupled Schwinger-Dyson equations for the scalar masses. Substituting 
the solution x((t) into V{a,x) gives the hnite temperature effective potential as a function of a 
alone: H(x((T),cr) = U{a). We have carried out numerical computations of U{a) for a range of 
coupling strengths and temperatures. 


In Fig. m we present numerical results for the choice of couplings Ai(M,.)/ 167 r 2 = A 2 / 167 r^ ~ 0.25, 
where we have taken ~ 4.5u/A^. This value for the loop expansion parameters is large, but 
still in a range that keeps the scalar masses, Eq. (^), somewhat below Attv/N, and the Landau 
pole above the momentum range relevant to the problem. For various temperatures, we plot 
N~'^U{a) + (7r2/45) for a a range over which U is real. The purely temperature-dependent 
second term merely shifts the curves onto the same scale, preserving their shapes and ordering. 


For a values somewhat smaller than those shown in the curves, below where the Nambu- 
Goldstone propagator mass x vanishes, U{a) develops a complex part just as in the 0{2N‘^) case. 
However, with A 2 7 ^ 0 the propagator masses in our approximation do not include all terms of 
higher order in A 2 at leading order in \/N. Goldstone’s theorem, applied at a minimum of the 
effective potential, refers not to such propagator masses calculated at lower order, but instead 


9 

















U+cT^4 



Figure 3: The effective potential N~‘^U{a) for several temperatures, shifted by (7r^/45) for 
clarity, with all dimensionful quantities measured in units of v/N. A first order phase transition 
mnst occur below the temperature Tinflect = 1.93u/A^ where the potential maximum and minimum 
coincide. The field strength discontinuity is at least 0.7 v/N here, with parameters Xi{Mj. = 
4:.54v/N) = 40 = A2. 


to second derivatives in symmetry directions. Thus y can vanish closer to the origin than the 
symmetry-breaking minimum. With the parameters of Fig. Tj is close to unity over the a 
range shown; we thus avoid infrared divergences at the minimum. 

The strongly first order character of the phase transition is apparent from the curves, in which 
the symmetry-breaking local minimum first appears at a field value fTmin ~ 0.7n/A7, at temperature 
^inflect just below 2v/N. With decreasing temperature the minimum evolves to larger field values 
until it becomes the usual T = 0 symmetry-breaking vacuum with cimin = v/N. We have checked 
that as A2 decreases, (Tmin(T'inflect) also decreases, with Tjnfiect varying only slightly. This behavior 
is consistent with the second order nature of the phase transition in the A2 = 0 limit. 

We can gain some qualitative insight into the origin of first order behavior in the A2 7^ 0 theory 
by contrasting it with the A2 = 0 theory. There, all but one of the degrees of freedom are Nambu- 
Goldstone modes, whose mass y vanishes at the potential minimum The vanishing of y 

leaves Eq. (|^ at iTmin in the simple form (7 ^j„ — {v/N/^ + T^/O = 0, clearly showing second order 
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behavior at T = Tc ^ ^/Qv/N. The relatively small masses generated by standard model gauge 
interactions would in fact turn these modes into pseudo-Nambu-Goldstone bosons, but such small 
perturbations induce at most weakly first order behavior. Even without those interactions, for 
small nonzero A2 the — 1) cja modes may be thought of as pseudo-Nambu-Goldstone bosons of 
the 0{2N‘^) theory, and their masses do not vanish with the VTa mass. For large A2, however, this 
is a large perturbation, preventing the equation for fTmin from taking on the simple second order 
form above. The full transcendental equation in our case leads to strongly first order behavior. 

4 Conclusions 

We have argued using an effective linear sigma model that the high temperature electroweak phase 
transition, in the presence of a strongly interacting symmetry-breaking sector, is strongly first 
order. With the global symmetry of this sector taken to be U{N) x U{N), the effective linear 
theory contains two dimensionless coupling constants Ai and A2, which are taken to approach the 
strong coupling limit. The range of couplings considered avoids the Landau pole problem, and 
corrections to the linear theory, represented by higher-dimensional operators, should be relatively 
small. Equivalently, a finite ultraviolet cutoff on the theory can be taken to he above the momentum 
range of interest. When the strong coupling limit is reached (loop expansion parameters equal 
unity), higher-dimensional operators are not obviously suppressed and the scalars described by the 
linear model are not obviously lighter than other new physics. Nevertheless, the strongly first order 
character of the transition as the couplings of the linear theory approach unit strength provides 
evidence that a theory such as technicolor will indeed give rise to a strongly first order electroweak 
phase transition. 

A 1/N expansion provides the framework for the computation. To leading order in 1/N, all 
contributions from the Ai interaction can be summed. With respect to A2, all planar graphs 
contribute to leading order in 1/N, and only a restricted set of contributions were summed. We 
argued that this is adequate to determine the strongly first order character of the phase transition. 
The strength of the transition decreases with A2, consistent with a second order transition in the 
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0{2N‘^) theory. 


The momentum scales in the strongly first order case are such that there is no reason to 
expect that a high temperature expansion should be reliable, or that the phase transition should 
be governed by an effective three-dimensional theory. Our results are nevertheless consistent with 


the hints provided by analysis of such theories |11, 151. 


This work may be relevant to the problem of electroweak baryogenesis. More detailed calcu¬ 
lations will be necessary to study the dynamics of departure from thermal equilibrium during the 
transition, but our strongly first order transition provides at least a necessary condition for baryon 
number generation. Furthermore, the sphaleron mass in our model is the same as the minimal 
standard model formula |^], msphaieron = K/g){ao), where iF ~ 3 is a function of the scalar 

couplings. Since is of order 2vjN, a strongly first order phase transition to a broken phase 
with {ao)/v ~ 1 ensures Boltzmann suppression of baryon number destruction during and after the 
phase transition. 
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